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Frustration introduces a nontrivial dispersion relation of spinwave even in a ferromagnetic phase in a spin system. 
We study the reflection and refraction process of spinwaves in the ferromagnet/frustrated ferromagnet junction system 
by using the Holstein-Primakoff spinwave expansion and taking the large-5 Umit. We discuss the relation between the 
incident angle and refraction angle of spinwave, namely, the Snell's low of spinwave. As concrete examples of frustrated 
ferromagnets, we study the fully polarized ferromagnet phases in the J1-J2 chains and the J1-J2 model on the square 
lattice. The interesting refraction processes, such as the splitting of the incident spinwave and the negative refraction, are 
discussed. We also study the transmittance and reflectance in these concrete models. 
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1. Introduction 

Spinwave is a collective excitation of spins in magnets. 
Since spinwave can carry information on the nano/microscale, 
recently, potential application of spinwave to new devices 
has attracted much attention to researchers. ^ In many circum- 
stances, spinwave behaves like a sound or light wave; reflec- 
tion and refraction process of spinwave can be considered. 

Understanding of the reflection and refraction process 
is important to control spinwave. Theoretically, this pro- 
cess of junction system of usual ferromagnets whose spin- 
wave has a trivial dispersion relation is understood by us- 
ing geometrical-optics approximation^"^ or Landau-Lifshitz- 
Gilbert equation.^ One of the important results of these stud- 
ies is that the direction of the propagating spinwave depends 
on the strength of applied magnetic field. The experimental 
control of propagation of spinwave is indeed realized by tun- 
ing the internal demagnetizing fields in a permalloy waveg- 
uid.^ If the dispersion relation of spinwave is nontrivial, an 
exotic refraction process is expected. In the case of light, for 
example, the negative refraction occurs due to the nontriv- 
ial dispersion relation.^' ^ Recently, it is theoretically proposed 
that the anisotropic nature of the dispersion relation of dipole- 
exchange spinwaves can result in the negative refraction.^ 

Recent experimental advance has been exhibiting a wide 
variety of frustrated magnets. Magnetic frustration provides 
us with an exotic magnetic behavior. Even in a ferromagnet, 
in which all spins align, the frustration can induce an unusual 
dispersion relation of spinwaves. Hence, it will be natural to 
expect interesting refraction processes of spinwaves in frus- 
trated ferromagnets. In this paper, we study the reflection and 
refraction process of spinwave in ferromagnet/frustrated fer- 
romagnet junction system. 

As a concrete example of frustrated ferromagnets, we 
consider the systems which are already realized experi- 
mentally. For example, frustrated ferromagnets can be pre- 
pared by applying high magnetic field on frustrated mag- 
nets. One-dimensional (ID) J1-J2 Heisenberg chain with 
a nearest-neighbor exchange coupling Ji and next-nearest- 
neighbor coupling J2 is one of the famous frustrated mag- 



the wavevectors Qs of cosQ = -Jil^Ji for I/1I//2 < 4 and 
J2 > 0. There are many compounds consisting of the J1-J2 
chains (see Table. I in Ref. 12). The saturated ferromagnetic 
phases under high magnetic field are experimentally realized 
ini2Rb2Cu2M030i2/^-^kN2H5)CuCl3 and22-24LiCuV04. 

As another example of the frustrated magnets, the var- 
ious compounds of the square-lattice J1-J2 model, e.g., 
BaCdVO(P04)2, are reported. ^^'^^ Without external field, the 
collinear antiferromagnet phase appears in BaCdVO(P04)2; 
the saturation field is about 4~6T. In the saturated ferromag- 
netic phase of this compound, theoretically, the dispersion 
relation is expected to be nontrivial; it has two minima at 
Q = (0, tt) and Q = (7r,0).25-27 

The organization of the present paper is as follows. In 
Sec. 2, we briefly review the dispersion relations of spinwaves 
in (frustrated) ferromagnets by using the Holstein-Primakofl' 
spiwave expansion. In Sec. 3, we generally discuss refraction 
angles of transmission spinwaves in the ferromagnets junc- 
tion systems satisfying the given conditions. Then, we apply 
this discussion to the several junction systems including frus- 
trated J1-J2 chains or J1-J2 model on the square lattice. We 
shall explicitly see the splitting of the incident spinwave or 
the negative refraction in these models. In Sec. 4, we study 
the reflection- and transmission rate (reflectance and transmit- 
tance) in some of the junction systems considered in Sec. 3 by 
solving the simple Schroedinger equations within the large S 
limit. We also study the phase shift of the spinwave at the 
boundary. Our approach naturally treats the lattice structure 
of ferromagnets without the coarse graining. 

2. Dispersion relation of Frustrated ferromagnets 

As a brief review, let us discuss the dispersion relation 
of spinwaves in the fully-polarized ferromagnetic phase of 
(frustrated) spin systems. For simplicity, we consider the lat- 
tice systems with one magnetic ion per unit cell, and assume 
the rotational symmetry around the z direction in spin space. 
This assumption leads to the conservation law of the total an- 
gular momentum along the z direction, namely, 2]i<5'5). We 
study the Hamiltonian with the generic exchange interactions 



nets.ii-15 In the fully polarized phase under high magnetic = ^i-J = ^J- = ^i' on-site anisotropic interaction term 



field, the dispersion of spinwave has two minima at 



11,16-18 
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and the external magnetic field H\ 

<i,j> i 



(1) 



where we use the coordinate {x,y,z) in the spin space, and 
(a, /?, c) in the lattice space. Let us rewrite the spin operator 
by using the Holstein-Primakoff' transformation on the fully- 
polarized ferromagnetic phase: 

^5 = -S +a]ai , 



(2) 



By using the 1 jS expansion, we approximately obtain the free 
bosonic Hamiltonian in the leading order in IIS: 




where 



7^ = ^(a;(k)-yu)a^Qfk, 

k 

m = ^ ^ /i cos k . 1 , ^(k) = 2S (6(k) - e^,) , 



(3) 



(4) 



= 2^(6(0) -6n,m)-(H + 2^ /r). 



and 6min is the minimum of 6(k). The conservation law of the 
angular momentum along the z direction leads to that of the 
total magnon number. For any exchange interactions, a suffi- 
ciently large H or ^ can induce the gap (ju < 0) in the magnon 
dispersion. If yu < 0, the ferromagnetic phase is stable (with 
few exceptions for finite 5"). Throughout this paper, we focus 
on this free bosonic Hamiltonian. 

We have seen that frustration can induce a nontrivial dis- 
persion relation. The nontrivial dispersion relation can lead to 
the nontrivial group velocity not parallel to the phase velocity, 
where the group velocity is given by 



v,(k) = Vk^(k) . 



(5) 



Physically, the group velocity carry the angular momentum as 
reviewed in Appendix. Not the phase- but the group velocity 
determines the traveling direction of spin wave. Next, let us 
study the various dispersion relations in the concrete models. 

2. 1 simple case 

First, we consider the simple Heisenberg model on the cu- 
bic lattice with the nearest neighbor ferromagnetic coupling 
/ < as shown in Fig. 1. 6(k) is given by 

6^(k) = J(cosk^ -h cos^^ -h cos^"") , 

(6) 

Wn = 6,(ko = (0,0,0)) = 3/. 

This leads to the dispersion relation ^^^(k) shown in Fig. 2. 
In the long wavelength limit |k| 0, this dispersion relation 
becomes isotropic and is given by 



k^ 

2m, 



1 
J 



(7) 



In the following discussion, we use this model as the concrete 
example of 'usual' ferromagnets. 




Fig. 1. (Color online) The cubic lattice with the nearest neighbor Heisen- 
berg exchange coupling /. The dots represent spins. The lattice constant 
flo = 1 is assumed. 




Fig. 2. (Color online) The dispersion relation ^^(k) of the trivial ferromag 
net on the cubic lattice for / < 0. = is assumed. 



2.2 J1-J2 chain 

Next, let us discuss the dispersion relation of the ID J1-J2 
chains with the ferromagnetic interchain coupling J3 < on 
the cubic lattice^^"^^'^^ as illustrated in Fig. 3: 

€i(q) = J I cos -h J2 cos 2q^ -h /3(cos q^ -h cos q^) 

= 2J2(cosq' -h —^f - 3^ - -^2 + Jsicosq"" -h cos^^) , 

(8) 

where the J1-J2 chains are assumed to lie parallel to the c axis. 
For I/1I//2 < 4, /2 > and J3 < 0, the dispersion relation has 
two minima at Q = (0, 0, ±Q^) where 

COSe^ = --f , 6nnn = -^-/2+2/3. (9) 
4j2 OJ2 

The dispersion relation 6i;i(k) is graphically shown in 
Figs. 4,5. 



Fig. 3. (Color online) Ji -J2 chains with the interchain coupling J3 on the 
cubic lattice. The dots represent spins. 
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Fig. 7. (Color online) The dispersion relation 6t>>2(k) for /1//2 = -0.5 and 
/2 > in the J\ -J2 model on the two-dimensional square lattice. 



1.0 




Fig. 5. (Color online) The dispersion relation cji (k) of the ID /i -J2 chains 
on the cubic lattice for /1//2 = 0-5, /2 > and J^IJ2 = -2. = is 
assumed. 



2.3 J1-J2 model on the square lattice 

The fully saturated phase in J1-J2 models on the square 
lattice (see Fig. 6) also has the nontrivial dispersion relation 
since 

e2{q) = Jiicosq"^ + cosq^) + J2(cos(q'^ + q^) + cosC^"" - q^)). 

(10) 

For -2 < /1//2 < 2 and J2 > 0, the dispersion relation 0J2(k) 
has two minima at Q^^^^ = (0, tt) and Q2^^ = (tt, 0) as shown in 
Fig. 7. 
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Fig. 6. (Color online) Two dimensional square lattice. The dots represent 
spins connected by the nearest neighbor Heisenberg exchange coupling J\ 
and the next nearest neighbor coupling J 2 . 



3. The Snell's low in the ferromagnets junction system 

In this section, we study the relation of the angles of 
incident-, reflected- and transmission- spin waves passing 
through a boundary (the Snell's law). Let us consider a ferro- 
magnet/ferromagnet junction system whose boundary plane 
is flat and is perpendicular to c direction. On the boundary 
the proximity efl'ect may induce a spin-exchange coupling be- 
tween spins on the boundary surface. This proximity efl'ect 
leads to transmission of incident spinwave. The Schematic 
figure of the reflection and refraction process in the usual fer- 
romagnets junction system is shown in Fig. 8. 











a 








' 



Fig. 8. (Color online) Schematic figure of refraction and reflection process 
of spinwave in the usual ferromagnets junction system in the case of k}) - 0. 
The travelling direction of spinwave is determined by the group velocity. 



In the case of light, the Snell's low is understood by the 
conservation low of the energy and that of the momentum 
parallel to the surface. The same context can be applied to 
determines the Snell's low of spinwave in the ferromagnets 
junction system if the following conditions are satisfied: (i) 
each ferromagnet and the exchange interactions on the bound- 
ary have a rotational symmetry around the z direction in spin 
space; (ii) each magnet have a translational symmetry; (iii) 
the lattice vectors in the a-h plane are same in both magnets. 
The assumption (i) leads to the conservation of a total num- 
ber of magnons, and (ii),(iii) lead to the conservation law of a 
total momentum of the a and h direction. Of course, the total 
momentum conservation law of the c direction does not hold 
by the boundary eff'ect. 

Since the above discussion is too formal, let us study the 
Snell's low of spinwave by considering a simple model. For 
example, one of the simplest model satisfying these condi- 
tions is the junction system consists of ferromagnets with the 
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nearest neighbor interaction / on the cubic lattice (see 
Sec. 2.1) with the same lattice constant. This setting is shown 
in Fig. 9. It is appropriate to consider the exchange interaction 
on the boundary as 

a,b 

where A is the strength of the exchange interaction on the 
boundary, and ci and C2 are the c components of the bound- 
ary's position. This boundary condition satisfies the condi- 
tions (i) and (iii). As long as (i) and (iii) holds, the detailed 
boundary condition aff'ect not the Snell's low of the spin wave 
but the transmission rate and the phase shift at the boundary. 
The explicit eff'ect of the boundary condition will be discussed 
in Sec. 4.1. 




Fig. 9. (Color online) Junction system of the usual ferromagnets on the 
cubic lattice. The dots represent spins. The exchange coupling A due to the 
quantum proximity effect is considered on the boundary. 



Let us consider the case that the incident spinwave in the 
medium (1) has the momentum k/. The incident spinwave ar- 
riving at the boundary is divided into the transmission wave 
with and the reflected wave k^. 

The assumption (i) leads to the energy conservation 

Q = OJ(i)(ki) - yU(i) = aJ(i)(kr) - yU(i) = aJ(2)(kt) - 11(2) . (12) 

As discussed in Appendix, the group velocity corresponds to 
the current of magnons. Hence, it is appropriate in our setting 
to assume that the group velocity of the transmission (reflec- 
tion) wave has the same (opposite) direction along the c axis 
as that of the incident wave. Namely, 

vf'(kdl\vf'(kd\ = -vl'^'(K)/\vl'^'(K)\ = vf ^(k,)/|vf ^(k,)| . 

(13) 

In addition, the momentum conservation law parallel to the 
boundary leads to 

=k^^=k^ ^ = k^^=k^ . (14) 

These equation determines the k^ and k^, which leads to the 
Snell's law. In general, the transimission- and reflected spin- 
waves determined by eqs. (12), (13), (14) are not necessarily 
single valued. 

Finally, let us discuss the general relation in the reflection 
and refraction process of spinwaves if the condition (i) (ii) (ii) 
are satisfied. If we label spinwaves as k(l>,k(2),k(3)..., this 



reflection and refraction process can be written as the ket 

K''+Z^i-^'^Zci<|>)|0)' (15) 

j ' j 

where |0) is vacuum state and Qf^^'^^(k) approaches the free 
boson in the media (1,2) with the momentum k away from the 
boundary. By considering the cuboid discussed in Appendix, 
we obtain the magnon number conservation low: 

\A\'vfHkd + \Bj\'vfHk'J') = X \Cfvf''(^^) ■ (16) 

j j 

The concrete expression of B and C will be discussed in 
Sec. 4.1 by explicitly considering the boundary condition. 

Since we have discussed the procedure to find the Snell's 
low, let us discuss the detailed properties of the Snell's low in 
the various ferromagnets junction system. 

3. 1 usual case 

First, let us further study the Snell's low of the junction 
system of the usual ferromagnets on the cubic lattice with the 
nearest neighbor exchange coupling / < (/' < 0) as shown 
in Fig. 9. The dispersion relation on each ferromagnet is re- 
spectively given by 

ojsiik) = 2^(/(cosr cosk^ cosk') - 3J) , (17a) 
0Js2(k) = 2S (/(cos r -h cos -h cos k') - 3/), (17b) 

where jii < and yU2 < 0. yUi,2 are given in eq. (4). In these 
usual ferromagnets, the group velocity v^, which carry the an- 
gular momentum, is given by 

\f\k) = Vk6(k) = -2^/(sinr, sink^ sink') , (18a) 

\f\k) = -2Sr{smk\ siny^^ sink') . (18b) 

The sign of the group velocity y^^^^^'^'^ is always the same 
as the sign of k^'^'' for -n < k^'^'' < n. If we consider the 
incident wave with wave vector k/, the reflected wave with k^ 
and the transmission wave with k^ are uniquely determined by 
eqs. (12), (13), (14). We obtain 

cosk', =j(J cos k^ -h (/ - r)(cosk^ cosJ^) - 3(J - f) 



(19) 

where k'l\k'\ = k'/\k'\ for -tt < k' < n. This implies that, 
the larger the chemical potential 112 is, the slower the group 
velocity of the transmitted wave becomes. Since [ii^ depends 
on the external magnetic field H, we explicitly see that the 
propagation of the spinwave can be controlled by the external 
magnetic field as discussed in Refs. 2-5. 

If we consider the long-wavelength limit |k| 0, the 
Snell's low in this system becomes simple. The group veloc- 
ity is proportional to the phase velocity in the long wavelength 
limit k ^ 0. Up to the order of k^, the refractive index is given 
by 

^sin^^^^ \ ^-r+li2 
^ sinOt \kc\ ^ Q-/-hyUi ' ^ ^ 

where 6i (6t) is the angle of incidence (refraction) as shown in 
Fig. 8, and Q is given in eq. (12). Eq. (20) is exactly the same 
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as the result obtained in Ref. 3. 

3.2 J1-J2 chains parallel to the boundary 

Next, let us consider the reflection and refraction process 
of spin waves in the usual ferromagnet/frustrated ferromagnet 
junction system. In this subsection, as a concrete example of 
frustrated ferromagnets, we consider the /1-/2 chains on the 
cubic lattice discussed in Sec. 2.2. We assume that the J1-J2 
chains are parallel to the boundary plane as shown in Fig. 10. 

The dispersion relation and the group velocity in the /1-/2 
chains on the cubic lattice are given by 

ojnO^) =2S(Ji cosr + /2Cos2r 



















' 



Fig. 11. (Color online) Schematic figure of refraction and reflection pro- 
cess of spinwave in the case of the negative refraction index and = 0. The 
travelling direction of spinwave is determined by the group velocity. 



.(11) 



(k): 



+ /(cosy^^ + cosyl" - 2) - - /2) , 

8/2 

: - 2^ (/i sin le + 2/2 sin 2^ , / sin / sin k') . 

(21) 



If we consider the incident wave of the wavevector k/ in the 
usual ferromagnet, the transmission wave of k^ is uniquely 
given by eqs. (12), (13), (14): 



cos^^ 



1 . ^+y^2 
J'^ IS 



J I cos^^ -/2COs2^;^ 



(22) 



-ncosl^^-\)--\--J2). 



where k^J\k^t\ = k^i/\k^\ for -tt < Z:^ < tt. As discussed in 
Sec 2.2, if -4 < /1/I/2I < 4 and J2 > 0, the equation 
Ji cos k^ + J2 cos 2k^ has two minima at cos k^^ = - ^ . If 

cosk^ > the a component of the group velocities of 

incident spinwave is of the opposite sign to that of the trans- 
mission spinwave: the refractive index becomes negative^^ as 
shown in Fig. 11. 







J 






1 1 

* i 


J 

i — i 


1 — ( 








1 — 1 








1 — ( 
1 — 1 












the /1-/2 chains is given by 

a;i2(k) = 2S (Ji cos k"^ + J2 cos Ik"^ 

A ^^^^ 

+ /(cosr + cos^^ - 2) + + J2) , 

8/2 

If the spinwave of the wavevector k/ in the usual ferromagnet 
is injected into the J1-J2 chains, the transmission waves of 
k^-^ are given by eqs. (12), (13), (14): 



cos^, 



.(±)c 



\Q.-2S /'(cos k^ + cos k^ -2)+ 112 Ji 

4/2 ' 

(24) 



45/2 



where = ut\ • This result implies that m 

some cases the two species of the transmission spinwaves ap- 
pear as shown in Fig. 13. 
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Fig. 12. (Color online) Junction system of the usual ferromagnet (left, 
green) /ferromagnet consisting of /1-/2 chains which lie along the c axis 
(right, red). The exchange coupling A due to the quantum proximity efl'ect is 
considered on the boundary. 



Fig. 10. (Color online) Junction system of the usual ferromagnet (left, 
green) /ferromagnet consisting of /1-/2 chains which lie along the a axis 
(right, red). The exchange coupling A due to the quantum proximity efl'ect is 
considered on the boundary. 



3.3 J1-J2 chains perpendicular to the boundary 

In this subsection, we study the case that the frustrated 
ferromagnet consists of J1-J2 chains is perpendicular to the 
boundary plane as show in Fig. 12. The dispersion relation in 



3.4 J1-J2 model on the square lattice 

Finally, let us briefly study the fully polarized phase in the 
J1-J2 model on the square lattice as an example of frustrated 
ferromagnets. We study the 2-dimensional junction system of 
the a-b plane as shown in Fig. 14. The dispersion relation in 
the J1-J2 model is given by eq. (10). In the case that -2 < 
J1/J2 < 2 and J2 > 0, by injecting the spinwave of k/ from 
the usual ferromagnet, the transmission spinwave of k^ in the 
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J 


' 



Fig. 13. (Color online) Schematic figure of refraction and reflection pro- 
cess of spinwave in the case of = 0. The incident spinwave can split into 
two transmission waves. The travelling direction of spinwave is determined 
by the group velocity. 



J1-J2 model is given by 

1 



COS 



/i +2/2 COS ^-^^ 2S 



Q + A/2 (25) 



where - -^i sin +2/2 cos sln^^ _ 
|/i sin qh+2J2 cos qa sin^/,| 



sinr 



^^7^. As well as the junction 

system discussed in Sec. 3.2, the negative refraction can be 
realized as shown in Fig. 11. 
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Fig. 14. (Color online) Junction system of the usual ferromagnet (left, 
green) /ferromagnet consisting of /1-/2 model on the square lattice (right, 
red). The exchange coupling A due to the quantum proximity eff'ect is con- 
sidered on the boundary. 



4. Reflectance and Transmittance 

We have studied the relation between the refraction angle 
of transmission spinwaves and the angle of incident spinwave, 
namely, the Snell's law of spinwaves in ferromagnets junc- 
tion system. In this section, let us discuss the reflectance and 
transmittance of spinwaves. Although in this section we fo- 
cus on the concrete models discussed in Sees. 3.1, 3.2, 3.3, 
our approach will be easily applied to the other frustrated fer- 
romagnets junction systems satisfying the conditions (i), (ii), 
(iii) discussed in Sec. 3. 

4. 1 usual case 

Let us study the reflectance and transmittance of spinwaves 
in the usual ferromagnet/ferromagnet junction system on the 
cubic lattice with the boundary condition (11) (see Fig. 9). As 
discussed in Sec. 3, the reflection and refraction process by 



the incident spinwave of k/ may be written as the ket 



a,b 



+cy.*'<-v,j|o), 



c=l 



(26) 



where we assume that the boundary is located between c = 
and 1, and is given by eq. (19). Then, A, B, C is determined 
by the Schroedinger equation: 



H\RRi} = Q\RRi) . 



(27) 



(28) 



(29) 



(30) 



(31) 



Obviously, (27) is satisfied at any lattice cite except at c = 
0, 1 . Hence, we consider 

{0\aa,b,c=oH\RRi} = ^(A + 5)^'W«+^f^) , 

{0\aa,t,c=iH\RRi} = ^lCe'^'y^'>^ . 
Equivalently, 

{eiikl) - Je-'^'i)A + {eiik:) - jS)B = AC . 
A(A + B) = (€2(k',) - rS)c , 

where 

el{k';) = 2J{co^k'.-\)-J + ^, 
€2(k',) = 2r(cosk', - 1) - / -h A . 
Hence, we obtain 

B _ (€i(k^) - Je-^^')(€2(k',) - re^^') - A^ 
A ~ (6i(y^p - Je^^i){e2m) - re'^') - A^ ' 
C liJ^mk": 
A ~ ~ (6i(y^p - Je^^'i)(€2(k^) - re'^t) - A^ ' 
By using the relation 

\{ei(l^) - Je-^'')(€2(k',) - rS) - A^p 
- 1(61 (i^p - jS){e2{K) - rS) - A^p 
= 4//A2 smk": sink', = Ahf\ki)vf\kt)/S^ , 

the magnon-number-conservationlaw given by eq. (16) is eas- 
ily confirmed as 

v':'HW-\j\') = vfHkr)\j\\ (33) 

regardless of any k/ and k^. From eq. (31), we explicitly find 
the dependence of the transmittance and the phase shift of 
spinwave on the boundary condition A. For A ^ 0, the trans- 
mittance leads |C/Ap = 0(A^/r^). 

4.2 J1-J2 chains parallel to the boundary 

Next, let us study the case that one side of the ferromagnets- 
junction system consists of the J1-J2 chains parallel to the 
boundary as shown in Fig. 10. The reflection and refraction 



(32) 
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process is described by the following ket: 

a,b c=-oo c=-oo 

oo 

C=l 

(34) 

where is given by eq. (22). The procedure to determine B, C 
is exactly the same as that of the previous section. Hence, B 
and C are given by (31). 

4.3 J1-J2 chains perpendicular to the boundary 

Finally, we consider the /1-/2 chains perpendicular to the 
boundary as a example of frustrated ferromagnets as show in 
Fig. 12. In this case, the ket is given by 

a,b c=-oo c=-oo 

(x> 00 

c=l c=l 

(35) 

where ^J*^"^ are given by eq. (24). We see that (27) is satisfied 
at any lattice cite except at c = 0, 1, 2. Hence, we consider 

{0\aa,b,c=oH\RRi} = Q(A + , 

{0\a,,b,c=iH\RR,} = Q(Ci + C2)e'<*°''^*'''> , (36) 

{0\aa,b,c=2H\RRi) = i2(Cie'*i''' + Cze'*!''')^'^*?''-'*?'" . 
Equivalently, 

(e[ik-:) - Je-"^)A + (6[(k'i) - Je"^)B = A(Ci + C2) , 
(6^(fc«^)-/i^*!^"-y2e'-2*!"')Ci 

+ ie^ik';'^') - Jie"^~" - /2e''*'"*")C2 = A(yl + B) , (37) 

+ ie'.ik',^^') - 2/1 cosk^-^ - /2e'''!"")C2e*"' = , 
where 

e'likf:) = IJicosk'j - I) - J + A 

72 

e^ikl^^') = e^(k^f >') = 2(71 cosA:,^^^" + J2 cos2A:,(+>" + -!-) 

0J2 

-J1+J2+A, 
e^(0-f2(0 + -^i-A. 

(38) 

B/A, Ci/A, C2/A are given by solving these equations. We nu- 
merically confirmed the following magnon number conserva- 
tion law in the various concrete parameters: 

)^(kO(l-|^P) = vf )^(kr^)|^^vf >^(k^>)|^|^ , (39) 

where v^^^^''(k) = -25 (/i sin^^ + J2 sin 2^^). 



5. Conclusion 

We studied the reflection and refraction process of spin- 
waves in the ferromagnet/frustratedferromagnet junction sys- 
tem whose ferromagnets are described by the spin Hamilto- 
nian with the generic Heisenberg-exchange coupling, the uni- 
axial anisotropic interaction and the external magnetic field. 
By using the Holstein-Primakoff' spinwave expansion and tak- 
ing the large S limit, the ferromagnet phase is given by 
the free bosonic Hamiltonian, which describes the dynam- 
ics of magnons. If frustration exists, the dispersion relation 
of magnons (spinwave) can become nontrivial. 

In Sec. 3, we discussed the Snell's low of the spinwave in 
the case that the following condition is satisfied in the junction 
system: (i) each ferromagnet and the exchange interactions on 
the boundary have a rotational symmetry around the z direc- 
tion in spin space; (ii) each magnet have a translational sym- 
metry; (iii) the lattice vectors in the boundary plane are same 
in both magnets. By studying the various junction systems 
satisfying these conditions, we found the nontrivial refraction 
process, e.g., the splitting of the spinwave and the negative 
refraction. 

In Sec. 4, we studied the 'reflectance', 'transmittance' and 
'phase shift at the boundary' in the concrete examples of the 
junction systems. Within the large S limit, we exactly ob- 
tain these quantities, which explicitly depend on the boundary 
condition. Throughout this paper our results does not need the 
long wavelength approximation. 

We thank G. Tatara, S. Murakami and A. Yamaguchi for 
useful discussions. 

Appendix: Current of angular momentum 

In this appendix, we briefly review the correspondence be- 
tween the group velocity of spinwave and the current of an- 
gular momentum. We study the spin system described by the 
Hamiltonian (1). Now, the Hamiltonian (1) commutes with 
the total spin S^^^ = Z/ Hence, S^^^ is the conserved quan- 
tity: the total magnon number 2i ^J^i conserves. This implies 
that the conserved current of magnons (related to the U(l) 
symmetry) exists. 

The magnon current is understood by studying in- and out- 
flow of the magnon number in a given region. Let us consider 
the cuboid Ci of which two face-to-face surfaces Si and 5 2 
are infinitely large as shown in Fig. A-1. 




Fig. A«l. Magnon current of the group velocity Vg passes through the 
cuboid Ci . 



The spatially-averaged flow of magnons is determined by 
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the current passing through S i and 52. By using the relation 



[4ai + a!aj,Si-Sj] = 0, 



(A-1) 



we obtain 

AZleC, rV t z^n c V r / t U 

I — = [2^alauH]=S 2j ^^1^^ - ai^-) , 

(A-2) 

where we use eq. (2). We assume that /y = Ji is short range so 
that J\ = for a large |1|, and that S\ is sufficiently separated 
from 52. Then, if we consider the state |k> = a^|0>, we obtain 

<k| '^^"aT Zj ^ysmk-(i-j) 



1 



= j^(As,nx • v^(k) + As,n2 • v^(k)) , 

(A-3) 

where A^. is the area of Si and 

v^(k) • n = 25 • n = -25 ^ /i sin(k • 1)1 • n . (A4) 

These equations imply that the group velocity can be viewed 
as the magnon current which carries the angular momentum. 
Even if there is a superposition of various spinwaves such 
as Z/A/a^ |0), the time-averaged magnon current is simply 
summed as 



(A-5) 
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